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The Polar Regions of Cassiopeia A: The Aftermath of a Gamma 

Ray Burst? 

J. Martin Laming^, Una Hwang^, Balint Radios^, Gergely Lekli^, & Endre Takacs^ 

ABSTRACT 

Probably not, but it is interesting nevertheless to investigate just how close 
Gas A might have come to generating such an event. Focusing on the northeast 
jet filaments, we analyze the polar regions of the recently acquired very deep 
1 Ms Ghandra X-ray observation. We infer that the so-called "jet" regions are 
indeed due to jets emanating from the explosion center, and not due to polar 



^ I cavities in the circumstellar medium at the time of explosion. We place limits on 

I the equivalent isotropic explosion energy in the polar regions (around 2.3 x 10^^ 

^ I ergs), and the opening angle of the x-ray emitting ejecta (around 7 degrees), 

O ! which give a total energy in the NE jet of order 10^° ergs; an order of magnitude 

^ I or more lower than inferred for "typical" GRBs. While the Gas A progenitor and 

"^H I explosion exhibit many of the features associated with GRB hosts, e.g. extensive 

presupernova mass loss and rotation, and jets associated with the explosion, we 
O I speculate that the recoil of the compact central object, with velocity 330 km 

^ \ s~^, may have rendered the jet unstable. In such cases the jet rapidly becomes 

^ ' baryon loaded, if not truncated altogether. Although unlikely to have produced 

a gamma ray burst, the jets in Gas A suggest that such outflows may be common 
^ ' features of core-collapse SNe. 

Subject headings: ISM: jets and outflows — supernova remnants — supernovae: 
individual (Gas A) — gamma rays: bursts 



^E. O Hulburt Center for Space Research, Naval Research Laboratory, Code 7674L, Washington DC 
20375 

lamingOnrl . navy . mil 

^NASA/GSFC Code 662, Greenbeh MD 20771, also Department of Physics and Astronomy, Johns Hop- 
kins University, 3400 Charles St, Bahimore MD 21218 

hwangSorfeo .gsfc .nasa.gov 



•^Institute for Experimental Physics, University of Debrecen, Bern ter 18/a, Hungary, H-4026 



- 2 - 



1. Introduction 

The speculation that gamma-ray bursts, or a subset thereof, might be connected with 
core-collapse supernovae was initiated by the coincidence of SN 1998bw with GRB 980425 
(Galama et al. 1998), and has been reinforced in recent years by the connections between SN 
2001ke and GRB 011121 (Garnavich et al. 2003), SN 2003dh and GRB 030329, (Hjorth et al. 
2003), and SN 20031w and GRB 031203 (Malesani et al. 2004). Spectroscopy of the afterglows 
of GRB 011211 (Reeves et al. 2002, 2003; Rutledge & Sako 2003; Butler et al. 2005), and of 
GRB 030227 (Watson et al. 2003) suggest the presence of highly charged ions of Mg, Si, S, 
Ar, and Ca, but not Fe, Co or Ni, (though Ni may be present in GRB 011211). Claims of 
the detection of spectral lines in the afterglows of GRB 001025A and GRB 010220 also exist 
(Watson et al. 2002), though both these gamma-ray bursts lack definitive SN connections. 
More recently though, Sako, Harrison & Rutledge (2005) have questioned the statistical 
significance of these GRB afterglow line identifications. It is against the background of these 
exciting developments that we turn our initial attention in analysis of the 1 Ms Chandra 
observation of the Cassiopeia A supernova remnant to the polar regions, with their striking 
"jet-fike" morphologies. We are motivated to examine in detail the spectra of ejecta knots 
in the so-called jet regions with a view to determining whether these regions are indeed due 
to an asymmetric explosion, and not just due to cavities in an asymmetric circumstellar 
medium, and if due to jets, whether or not we can infer energetics and other parameters 
connected with the jet nature. 

Ideas that core-collapse supernovae, particularly those occurring after extensive presu- 
pernova mass loss as Type Ib/c, might be inherently asymmetrical have also gathered force 
following SN 1987A (see e.g. Wang et al. 2002, and references therein). Linear polarizations 
of optical light of order a few percent are typically seen, suggesting nonspherical scattering 
stellar envelopes, with aspect ratios as large as 2. Mechanisms by which core collapse su- 
pernovae become asymmetrical generally derive from the rotation of the progenitor, which 
as it collapses may or may not produce strong magnetic field by a magneto-rotational insta- 
bility. A comprehensive review of recent work in this area is given by Wheeler & Akiyama 
(2004). Such ideas are appealing in the context of Cassiopeia A since rotating core-collapse 
and the associated magnetic fields (e.g. Proga et al. 2003; Proga 2005; Akiyama et al. 2003) 
or anisotropic neutrino emission (e.g. Yamasaki & Yamada 2005) are often suggested as a 
means of producing bipolar outflows or jets. Before proceeding further in this direction we 
need to be sure that the morphology we refer to as "jets" really are due to a feature of 
the explosion and are not arising as a result of cavities in the circumstellar medium at the 
appropriate locations. Blondin, Lundqvist, & Chevalier (1996) model the latter situation, 
and elongated structures of supernova ejecta resembling jets are easily produced. Gas A is 
inferred to have undergone extensive mass loss from its original 20-25 Mq progenitor, to 
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have been only 3-4 Mq upon explosion (Laming & Hwang 2003; Chevalier & Oishi 2003). 
Additionally, the surrounding remnant stellar wind is relatively dense and slow moving. Such 
conditions do not generally arise with a radiatively driven wind from a 20-25 Mq progenitor 
(Woosley, Langer, & Weaver 1993), but require the existence of a binary companion to aid 
the mass loss. In such a case one might expect a departure from spherical symmetry in the 
wind, from the binary orbital plane to the rotation axis, as is apparent in the case of SN 
1987A (Sonneborn et al. 1998). 

In this paper we use spectra from the newly acquired deep Chandra observation to 
investigate the properties of the polar regions. Wc describe the data and methods of analysis 
in section 2. In section 3 we demonstrate that models of the polar regions based on a 
circumstellar cavity to explain the morphology cannot reproduce the observed spectra, and 
that we really are seeing the evolution into the remnant phase of a bipolar explosion. Section 
4 discusses the nature of these jets in more detail and section 5 concludes. 

2. Observations and Data Analysis 

We use spectra of the polar regions of Cas A from the Chandra Observatory Very Large 
Project 1 Ms observation. These data were taken with the backside-illuminated S3 CCD 
chip of the Advanced CCD Imaging Spectrometer (ACIS) in 9 observation segments, mostly 
over two weeks in 2004, as described in more detail in Hwang et al. (2004). The data were 
corrected for the time- dependent gain across the S3 chip, but could not be corrected for 
charge-transfer inefficiency because the photon events were graded on board the spacecraft 
prior to data telemetry. 

Several regions along the three main filaments of the NE jet were chosen for spectral 
analysis, as shown in Figure 1, as well as the composite of the two main filaments of the much 
fainter counterjet in the southwest (see Hwang et al. 2004). The southernmost filament in 
the NE is the straightest and longest of the three, and extends farthest into the interior of 
the remnant. Of particular interest are the very faint knots at the very outer tip of the NE 
jet filaments that are newly revealed by the deep Chandra observation. Their spectra are 
shown in Figure 2. 

All the spectra were extracted individually for each of the nine observation segments, and 
corresponding detector response files generated using CIAO v3.0.2 (http:www.cxc.harvard.edu/ciao). 
The individual spectra were then added together to produce the final spectrum for a par- 
ticular region, while the individual response files were weighted according to the relative 
exposure time of each observation segment before being added to produce the final response 
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files. The time-dependent accumulation of soft X-ray absorbing contaminant on the ACIS 
detector was modelled during spectral fitting with the ACISABS model component. The 
detector gain was adjusted to optimize the fits, compensating in part for uncertainties in the 
energy scale, due both to detector performance and the intrinsic bulk motion of the gas. 

The goal of the spectral analysis is to obtain electron temperature, ionization age (the 
integral of electron density over time since shock passage for the reverse shocked ejecta, 
assuming a constant electron temperature during this time), and element abundances. We 
follow Laming & Hwang (2003) in fitting the spectra with models for single-temperature 
plasmas with timc-dcpcndcnt ionization (nonequilibrium ionization or NEI). We also exper- 
imented with collisional ionization models, NEI models with multiple components, and NEI 
models with a range of ionization ages (mimicking plane-parallel shocks). For the knots 
along the main body of the jet filaments, the spectral fits were generally improved by using a 
plane-parallel shock (pshock) model compared to a single NEI component. Two-component 
NEI models were also an improvement, but generally did not perform better than the pshock 
models. It is noteworthy that the knots at the very tip of each jet filament have spectra 
that are well-described with a single component (although this is at least partly because 
these spectra contain fewer counts). The ionization ages fitted for the jet tip spectra are 
very high, and they may also be successfully fitted with collisional ionization equilibrium 
models (see Table 1). All the jet region spectra show strong emission lines of Si, and S. 
Those with sufficient counts also show Ar and Ca as well as Fe K emission. Except for the 
rather extended knot at the tip of the middle NE filament, the jet tip spectra are effectively 
cut off at around 5 keV. 

There are various possibilities as to which elements contribute to the continuum. Fol- 
lowing Laming & Hwang (2003), we generally model the continuum as coming from either 
ionized O or Si, with an O continuum consistently giving a better fit (those results are given 
in Tables 1 and 3). On the basis of alone, however, it is not possible to confidently 
distinguish between the various continuum models. The assumption of which light elements 
provide the continuum does not strongly affect the quality of the spectral fits but does affect 
inferences concerning the electron densities and masses in the knot. H continuum models, 
for example, require higher densities and masses than those with O continuum. 

Given the complexity of the spectral models required for most of the jet knots, it is 
helpful to examine properties of the line emission to assess qualitative trends in the data. 
In particular we are interested in possible trends that might be present with location along 
the jet in the line ratios or line equivalent widths. We therefore fitted the spectra with line 
blends of individual elements and a bremsstrahlung continuum, as summarized in Table 2. 
We illustrate in Figure 3 the two strongest trends that we found. The first panel shows 



- 5 - 



Fe K EQW plotted against the projected distance of the jet region from the point source 
near the center of the remnant. With the exception of the knots at the very tip of the 
jet (where the temperatures are significantly lower than elsewhere along the jet, limiting 
both the cmissivity of and sensitivity to Fe K), the Fe K line strength is seen to either stay 
roughly level or increase with distance outward along the jet filament. The northern filament 
in particular appears to show the strongest tendency for an outward increase in Fe K line 
strength. Although the line strength of Fe K is affected by a number of factors, this trend 
appears to be echoed in the fitted Fe abundances as well. The second panel of the figure 
shows the correlation between Ca He a. EQW and Fe K EQW. The Ca He a. blends in the 
jet knots are unusually prominent compared to other regions of the remnant (Hwang et al. 
2000), and are seen to roughly follow the Fe K line strength along the jet filaments. 

In the jet models discussed below, a significant fraction of reverse shocked jet plasma 
has now cooled by radiation and adiabatic expansion to temperatures such that it would 
no longer emit X-rays. The extrapolation of the observed equivalent width trends to larger 
distances outward along the jet is consistent with such models. If the element composition 
of the presently unseen jet tip material were dominated by Fe, as would be consistent with 
its origin deep inside the progenitor, then the radiative cooling time of such plasma would 
be even shorter than for the O dominated composition assumed so far, allowing it to cool 
very quickly. 

3. Circumstellar Cavities or Ejecta Jets? 

3.1. Circumstellar Cavities 

The bi-polar morphology exhibited in Cas A has two plausible origins. Both a symmet- 
ric explosion into an asymmetric circumstellar medium, or an asymmetric explosion into a 
symmetric medium could produce jet-like structures. The first possibility has been consid- 
ered by Blondin, Lundqvist, & Chevalier (1996). In core-collapse supernovae, asymmetries 
in the circumstellar density may arise from non spherically symmetric pre-supernova mass 
loss, as has been established in the case of SN 1987A. The high degree of mass loss in Cas 
A inferred from the dynamics (Laming & Hwang 2003), and the likelihood that such winds 
were driven by the interaction of the progenitor with a binary companion lend support to 
this idea. 

We simulate the evolution of the ionization balance for a knot with composition resem- 
bling that in the middle jet tip knot discussed above, i.e 0:Si:Fe being 0.82:0.13:0.05 by mass. 
The departure from spherical symmetry means that forward shocked circumstellar plasma 
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at the head of the jet moves non-radially (i.e. tangentially) away from the apex, with the 
consequence that the contact discontinuity at the apex is now much closer to the forward 
shock. Hence the forward shock in the polar region should be close to the outermost X-ray 
emitting ejecta, at about 3.8 pc, giving an aspect ratio of 3.8/2.5 = 1.5, taking the equatorial 
forward shock radius as 2.5 pc. Such an aspect ratio agrees best with a density contrast of a 
factor of eight between equatorial and polar regions (see Figs 3 and 4 in Blondin, Lundqvist, 
& Chevalier 1996). Accordingly we take our equatorial core-envelope model from Laming 
& Hwang (2003), with an outer envelope power law density profile with an exponent of 9 
and a uniform density inner core, and reduce the circumstellar density by a factor of 8 in 
our initial study. By following the evolution of the ionization balance and electron and ion 
temperatures in the ejecta following reverse shock passage, we compute, and plot in Figure 4, 
the locus of electron temperature against ionization age for this model, as well as for models 
with outer envelope power law exponents of n = 7, and 11. The upper branch of each curve 
corresponds to the phase when the reverse shock is propagating through the constant den- 
sity inner ejecta core. The lower branch corresponds to reverse shock propagation through 
the outer power law envelope, and the maximum ionization age is found at the uppermost 
point in each curve at the core-envelope boundary. For n — 9, this is found to be 4.3 x 10^^ 
cm~^s (log net = 11-6), with a temperature of 6.5 x 10^ K (logTg = 6.8). The other models 
either give higher ionization age and lower temperature, or the reverse. Comparison with 
the fit results from section 2, plotted as boxes corresponding to the uncertainties in temper- 
ature and ionization age for each knot, shows that no cavity model gives sufficient density 
of plasma at high enough temperature to match the observations. Simply put, the ejecta 
expand sufficiently rapidly into the cavity that the density is too low either for appreciable 
electron-ion equilibration to raise the electron temperature or to ionize the plasma to the 
observed values. The only way to make such a model work would be to include some col- 
hsionless electron heating at the reverse shock (currently neglected, see Laming & Hwang 
2003, and and references therein for fuller discussion of this point), which with reference to 
Figure 4 could provide the extra factor of three or so in electron temperature required to 
match the observations of knots in the jet "stem" , but those at the jet tip would still be 
discrepant. 



3.2. Jets 

We model jets based on the simulation of Khokhlov et al. (1999). These authors studied 
the explosion of the inner 4.1 Mq of a 15 Mq progenitor (assuming substantial presupernova 
mass loss), induced by baryonic jets emanating from the central regions of the star, where 
the composition is dominated by Fe and Si. The inner 1.6 Mq is assumed to fall back onto a 
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protoneutron star, and the remaining 2.5 Mq of ejecta should give a reasonable match to the 
parameters of Cas A. Initially, each jet comprises about 0.05 Mq of material and together 
they have 9 x lO^'^ ergs of kinetic energy. The jet power is constant for the first 0.5 s, and 
ramped down to zero after another 1.5 s. The equivalent isotropic mass and energy in the jet 
when they are launched would be 2 Mq and 2 x 10^^ ergs respectively, increasing to 8 Mq and 
8 X 10^^ ergs upon jet breakout. From the radius of the jet at launch compared to its final 
radius upon emergence through the stellar surface, we infer a jet density profile pj oc 1/r, 
which also matches the figures presented by Khokhlov et al. (1999). We use an adaptation of 
the BLASPHEMER code (Laming & Grun 2002; Laming & Hwang 2003) which implements 
the solutions for ejecta profiles p oc r~"^ with n <3 given by Truelove & McKee (1999), and 
is summarized here in Appendix A. A jet unconfined by the surrounding ejecta would obey 
Pj oc 1/r^. 

In Figure 5 we plot the locus of electron temperature and ionization age for jet models 
assuming 2.3 x 10^^ equivalent isotropic energy and 1.815 Mq ejecta mass. The ejecta mass 
is taken from that determined in the polar region exterior to the jet by Laming & Hwang 
(2003). The energy is modified slightly from the value used by Khokhlov et al. (1999) to 
give better agreement with the observed morphology. Specifically, we want the faint X-ray 
knots observed at the very tip of the jet to end up at a radius of approximately 3.8 pc from 
the explosion center. The models differ in that abundance sets of 0:Si:Fe are taken from 
the fits to the three jet tip knots, color coded as black - north tip, 0.004:0.628:0.368; red - 
middle tip, 0.82:0.13:0.05; and green - south tip, 0.55:0.12:0.33. Additionally we have also 
computed a model appropriate for the north tip, where half the plasma is assumed to be 
H. This is given by the black dotted line. As can be seen, the jet models give a reasonable 
match to both the CIE knots at the jet tip, as well as the NEI knots further down the 
jet stem, in terms of ionization age and electron temperature. Some of the jet stem knots 
fall at lower temperatures than modeled assuming O to be the most abundant clement, but 
which could be understood if the jet knot plasma comprises a substantial fraction of H. Of 
course the abundances in the jet stem knots are probably not the same as those in the three 
tip knots, although here they have been modeled as such. The variation between the three 
curves for N, M, and S jet filaments should give some idea as to the variation to be expected 
as abundances vary within reasonable limits. 

As well as comparing electron temperatures and ionization ages coming from our mod- 
els with those from fits, we have also investigated the ionization balance predicted by the 
models for Fe and Si in the jet tip knots. Figure 6 shows the model Fe and Si charge state 
distributions for ejecta that went through the reverse shock at various times between 1.3 
and 1.6 years after explosion, compared with the charge state distribution of a plasma in 
coUisional ionization equilibrium at 8.6 x 10^ K (the thin solid line), meant to match the 
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conditions in the jet tip knot of the middle filament. For both Fe and Si, the model charge 
states corresponding to an electron temperature of 8.6 x 10^ K are a poor match to the 
observed charge states, with too many ions in high charge states, e.g. Fe^^"*" and Si^^^. This 
is because the recombination from these ions is a slow process. Much better correspondence 
between observation and the models is found for ejecta that have been allowed to cool down 
to temperatures 5 — 6 x 10^ K, i.e. lower than actually observed. We speculate that although 
we have tried to model and fit the knot spectra as though they are single plasmas, shocked 
instantaneously by the reverse shock, some inhomogeneity may still exist. Small variations in 
elemental composition or density are likely to reveal themselves during the onset of thermal 
instability, when the temperature is decreasing very quickly with time, making the heavy 
element line emitting regions cooler than the surroundings which mainly emit continuum. 
At 1-2 years after explosion, the reverse shock speed in our models is 4000-5000 km s~^, 
which means it would take 10^ seconds or 3 years to traverse a knot of 1 arcsecond extent. It 
is therefore also likely that our assumption that the whole knot is shocked instantaneously 
breaks down, in that the reverse shock speed evolves as it traverses the knot making a single 
ionization age less realistic. 

A related problem is that knots in the main body of the jet are better fit by plane-parallel 
shock models that integrate over a range of ionization age, rather than by single ionization 
age models. These knots with smaller ionization ages are shocked later (i.e. around 5 - 
10 years after explosion) than those at the jet tip, at a time when the reverse shock in our 
models is evolving much less rapidly with time. These knots are also observed at smaller radii 
relative to those in the jet tip than our models would predict, and have similar problems with 
the ratios of He-like to H-like Si, in that the models predict much more H-like than actually 
observed. This is quite different from Si in knots of similar temperature and ionization age 
studied in Laming & Hwang (2003). The difference lies in the fact that for the jet models 
considered here, the electron temperature was previously much higher than it is at present, 
and the plasma is recombining. The NEI models which assume a constant temperature are 
probably not good approximations. This is quite different to the knots in more equatorial 
regions in Laming & Hwang (2003), where the postshock knot electron temperature is roughly 
constant, making the NEI models excellent approximations. We suspect that these knots may 
have initially undergone an oblique interaction with the reverse shock, leading to a lower peak 
electron temperatures, and have undergone further interactions with shocks refiected from 
the blast wave as it encounters clumps (i.e. quasi-stationary fiocculi) in the circumstellar 
medium, which would further decelerate and heat them. Such a scenario does indeed seem 
to be the case for SN 1987A (Zhekov et al. 2005), in that a variety of shock conditions 
are required to fit the Chandra/LETG spectra, and indeed should be expected as the blast 
wave begins to encounter the inner ring of circumstellar material, producing a succession of 
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reflected shocks back into the following ejecta. 

3.3. Related Issues 

Our conclusion that we are indeed seeing an asymmetric explosion into a symmetric 
medium is based solely on the spectroscopy of knots in the polar regions. Laming & Hwang 
(2003) using similar techniques for ejecta knots observed elsewhere in the remnant concluded 
that the explosion energy per unit solid angle was higher in regions close to the jet than in 
more equatorial directions. A completely independent conclusion of similar anisotropy in 
the explosion comes from a survey of the outlying optical knots by Fesen (2001). These 
optical knots are generally taken to be essentially undecelerated by their interaction with 
the circumstellar plasma, and hence any velocity variation exhibited by these knots must 
have an origin in the explosion itself, so that our conclusion is in line with those of other 
researchers. 

However considerable uncertainty surrounds the mcchanism(s) of formation of the vari- 
ous kinds of optical knots, so here wc briefly offer some speculations as to how these features 
fit in with the hydrodynamic models for Cas A studied in Laming & Hwang (2003). Knots 
rich in nitrogen are observed around the remnant (but not apparently in the NE jet region) 
expanding with velocity about 10,000 km s"^. Due to their composition, they must originate 
in the outer layers of ejecta. We suggest that these knots are the end result of a Rayleigh- 
Taylor instability at the contact discontinuity. In simulations, R-T fingers of ejecta typically 
do not penetrate through the forward shock, but Jun, Jones & Norman (1996) show that 
in the presence of a clumpy external medium, this is not the case. The R-T fingers can 
gain extra kinetic energy from the vortices generated by the shock-cloud interactions. In 
this case we would expect the velocity with which the resulting clumps of nitrogen rich ma- 
terial expand to be similar to the ejecta expansion velocity at the core-envelope boundary, 
which in the models of Laming & Hwang (2003) is indeed around 10,000 km s~^. Knots rich 
in oxygen and its burning products (they also show prominent S II emission) are found in 
various regions of the remnant, especially in the NE jet. Many authors have surmised that 
these arise as a result of the "Fc bubble" effect, whereby Ni-Co-Fc bubbles inflated by the 
energy deposition due to the radioactive decay of Ni to Co and ultimately to Fe compress 
the surrounding plasma. This argument would suggest that Fe should be present, or should 
have been present in the jet regions in order to produce the fast moving knots observed there 
today. We will return to the subject of Fe in the jets below. 

Finally, the location and velocity of the reverse shock in the NW and SW regions of 
the remnant has been inferred from optical observations with HST WFPC2 images taken in 
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2000 and 2002 (Morse et al. 2004). The reverse shock radius of ~ 1.8 pc and velocity with 
respect to the freely expanding ejecta of ~ 2000 km correspond reasonably well with the 
n — 9 model from Laming & Hwang (2003), which predicts a reverse shock radius of 1.7 pc 
and velocity of 1500 km s^^. The ejecta free expansion velocity in this model is 4800 km s~^, 
to be compared with 5000 km s~^ in Morse et al. (2004). Adjusting the model of Laming 
& Hwang (2003) to better match say the faster reverse shock velocity (e.g. by decreasing 
the ejecta envelope power law index to n — 7) would make the discrepancy between the 
measured and modeled reverse shock radii worse, and vice versa. We would expect the 
models to slightly underestimate the reverse shock velocity at 340 years since Laming & 
Hwang (2003) make the approximation of holding the reverse shock velocity constant during 
its propagation through the uniform density ejecta core, whereas in fact we would expect 
it to slowly accelerate. Consequently the reverse shock velocities computed at earlier times 
in the remnant evolution for a given model should be more accurate than those for the 
current epoch. In any case, the reverse shock is observed to be fragmented and irregular in 
morphology due to local density inhomogeneities, and may even be dynamically unstable, so 
we consider the agreement between our model and the HST observations to be satisfactory. 



4. Discussion 

4.1. Jet, Mass, Energy and Opening Angle 

We first discuss in more detail the nature of the Truelove & McKee (1999) n—1 models 
employed to simulate the jet emission. The equivalent isotropic energy of 2.3 x 10^^ ergs is 
chosen to place ejecta cooling through temperatures of 0.7 keV at a distance of about 3.8 pc 
from the explosion center. These models predict a blast wave radius at the jet tip of 5.66 
pc and a velocity of 7850 km s^^. For reasons to be discussed below, these are probably 
overestimates, but we do not actually detect the location of the blast wave anywhere in the 
jet region. This radius however is not too far off from the radius of the outermost optical 
knot in the jet of 4.8 pc (Fesen 2001). The gas that we see today at 0.7 keV temperature 
encountered the reverse shock between one and two years after explosion and has a mass 
coordinate of about 0.5. This means that nearly half the ejecta went through the reverse 
shock during these first one to two years, and has by now cooled to temperatures below those 
where it would be visible in X-rays to Chandra. Consequently in the jet there should be a 
lot of cold plasma at the jet head, cooled by both radiative losses and adiabatic expansion. 

To our knowledge, there are no published simulations of jet explosions carried through 
to the supernova remnant phase of the evolution. However we can get some idea of the 
hydrodynamics involved in the interaction of a jet with an ambient medium from recent 
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laboratory work (Blue et al. 2005; Foster et al. 2005; Lebedev et al. 2005). Here we can clearly 

see a jet stem, such as we plausibly observe in X-rays in Cas A, and a denser "mushroom 
cap" jet tip, which in Cas A has presumably cooled to too low a temperature to be visible in 
X-rays. From Figure 5 it appears that the knots at the jet tip are on the verge of becoming 
radiatively unstable. Radiatively cooled cjccta at the jet tip will reduce the pressure driving 
the forward shock, so that the forward shock radius and velocity might be expected to be 
smaller than in the adiabatic models of Truelove & McKee (1999). 

As mentioned above, we do not detect the blast wave in the NE jet region. Elsewhere in 
the remnant it is easily identifiable as a rim of continuum emission, presumably synchrotron 
emission from cosmic ray electrons. The width of this rim is identified as being due to the 
radiative loss time of the cosmic ray electrons as they are advected away from the shock 
(Vink & Laming 2003), with the result that the magnetic field must be of order 0.1 mG, 
presumably amplified from a very low value in the presupernova stellar wind by a cosmic 
ray precursor as suggested by Bell & Lucek (2001) and Lucek & Bell (2001). Why is such 
emission absent from the jet region? One possibility is that the shock here is oblique; the 
inflowing and outgoing plasma do not flow along the shock normal, and diffusive shock 
acceleration is known to be less efficient in such cases (e.g Drury 1983), leading to lower 
cosmic ray energies and weaker magnetic field amplification by the cosmic ray precursor. It 
is also possible that the blast wave at the jet tip is actually outside the Chandra field of 
view. 

Our simulations of the knot spectra above give the equivalent isotropic energy in the jet 
as 2.3 X 10^^ ergs, with the assumption that the jet equivalent isotropic mass is 1.815 M©, as 
determined from knots at the jet base by Laming & Hwang (2003). We are however, unable 
to independently constrain the jet mass in the same way since as mentioned above, we do not 
detect the blast wave driven by the jet. The plot of electron temperature against ionization 
age is essentially independent of the jet mass and energy, so long as these are constrained to 
place plasma cooling through a temperature of about 700 eV at a radial distance from the 
explosion center of about 3.8 pc. A higher jet mass (as in e.g. Khokhlov et al. 1999) would 
of course require a higher jet energy, but we consider this an unlikely possibility. We have 
made estimates of the plasma density and mass in the knots at the jet tip from the spectral 
fits, which generally come out lower than our model values. While there is considerable 
uncertainty in these estimates, they suggest that if anything the jet has density similar to 
or less than the surrounding stellar envelope, and is not overdense as in the case modeled by 
Khokhlov et al. (1999). 

Given an equivalent isotropic energy in the jet of 2.3 x 10^^ ergs, the evaluation of the 
total jet energy requires some knowledge of the opening angle. In Figure 7 we overlay blast 
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wave profiles including the NE jet estimated from a "pseudo" Kompaneets approximation 
outlined in Appendix B. We assume that sufficient cold ejecta, either shocked or unshocked, 
exists in the jet region so that the remnant retains some "memory" of the initial energy 
distribution. We take the jet energy profile to be constant within the jet opening angle, and 
to fall off as oc 9^^ outside this region. We consider values of n of 2, 3, and 8, coming from 
values favored by Lazzati & Begelman (2005), Zhang, Woosley, & Heger (2004) and Graziani, 
Lamb, & Donaghy (2005) respectively. The panels show the cases of jet opening angles of 
5 degrees (top left), 7 degrees (top right), 9 degrees (bottom left) and 11 degrees (bottom 
right), and in each case the n = 8 model is the narrowest and n — 2 the widest. The total jet 
energy evaluates to about 10^° ergs, taking an opening angle of 7 degrees. This is nearly an 
order of magnitude lower than the "standard" gamma-ray burst radiated energy found by 
Frail et al. (2001), and probably much lower than the kinetic energy of rclativistic cjccta in 
such events, and argues against Gas A actually having produced a gamma-ray burst during 
its explosion. This energy, with the likely jet underdensity compared to the surrounding 
stellar envelope also suggest that the jets are unlikely to have induced the explosion, being 
significantly less energetic than those in the simulation of Khokhlov et al. (1999). More likely 
they are the by-product of an explosion proceeding by different means. 

Our estimation of blast wave profiles around the jet region assumed that at all times in 
the evolution of Gas A there has been sufficient cold ejecta at the jet head to retain some 
"memory" of the initial explosive energy distribution with angle. In section 2 we described 
how the equivalent width of the Fe K feature is seen to increase with radial distance, especially 
in the northern filament, a trend which if extrapolated would suggest that the jet tip regions 
are dominated in composition by Fe. This would shorten the radiative cooling time for 
shocked ejecta over that calculated assuming O dominated composition, and support our 
assumptions here. Hwang & Laming (2003) estimate that only a few per cent of the total 
mass of Fc that Gas A is expected to have ejected arc currently visible as X-ray emitting 
ejecta, and so a significant portion of the unseen jet material could also be composed of Fe. 
A number of authors (Wang et al. 2002; Nagataki et al. 1998, 2003) have also suggested that 
bipolar jets should be locations where a-rich freeze out occurs with correspondingly large 
abundances of ^^Ti. In the X-ray emitting portion of the NE jet, we see rather little Fe, and 
certainly no knots that may be considered as pure Fe as found by Hwang & Laming (2003). 
We have searched for inner shell line emission from ^'^Sc and ^^Ga (the decay products of 
^"^Ti, formed mainly by K shell electron capture) in the region at the head of the jet where 
we expect the X-ray dark ejecta to be, without success. Our (3 sigma) upper limit is about 
4 photons, less than 0.1% of the "^^Sc inner shell line emission expected to be observed in 
the whole remnant, based on the gamma-ray line observations of nuclear deexcitation in 
^^Sc (Vink et al. 2001; Vink & Laming 2003) and ^^Ga (lyudin et al. 1994; Schonfelder et 
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al. 2000).-'^ We should caution that the regions we investigated are right at the edge of the 
Chandra/ ACTS field of view, and that regions that underwent a-rich freeze out may even be 
further out. 



4.2. Could Cas A Have Produced a GRB? 

As mentioned above, although Cas A was produced in an explosion that accelerated jets 
of material, the effects of which persist to this day, the energetics of these features make it 
unlikely that Cas A actually produced a gamma-ray burst. In this section we discuss how 
other features of the Cas A SNR might reinforce this conclusion. 

The Cas A explosion energy, in the range 2— 4 x 10^^, ergs places it clearly as a supernova 
and not a hypernova, though an observer placed on the jet axis (equivalent isotropic energy 
of 2.3 X 10^^ ergs) might have perceived it as such from line broadening (Mazzali et al. 2002; 
Ramirez-Ruiz & Madau 2004). Wang (1999) identified possible hypernova remnants in MlOl, 
supporting the suggestion (Paczynski 1998; MacFayden & Woosley 1999) that such a class of 
highly energetic explosions exists. Further, Atoyan, Buckley & Krawczynski (2005) identify 
a source detected only in TeV gamma-rays by HESS (High Energy Stereoscopic System) as 
the remnant of a gamma-ray burst, based on theoretical arguments. As Ramirez-Ruiz & 
Madau (2004) discuss, the absence of a Compton scattered jet component (i.e. scattering 
photons from a misaligned jet into an observer's line of sight) in otherwise normal Type 
Ibc SNe would suggest that those hosting GRBs are not normal Type Ibc SNe, supporting 
the conjecture that Type Ibc SNe hosting gamma ray bursts constitute a separate class of 
more energetic SNe, known as hypernovae. However more recently, Soderberg et al. (2006) 
present evidence that the optical properties of Type Ibc SNe associated with GRBs arc not 
significantly different to "normal" Type Ibc events, and that GRB associated SNe can only 
be rehably identified from their radio signatures. Nagataki, Mizuta & Sato (2006) show that 
^^Ni production in a GRB jet would be insufficient to power a hypernova, and that if these 
explosions exist as a separate class, ^^Ni must be produced elsewhere in the ejecta, further 
loosening any association between hypernovae and jets. 

Podsiadlowski et al. (2004) and van Putten (2004) estimate the branching ratio from 
Type Ibc supernovae to gamma-ray bursts to be in the range 10~^ — 10~^. Podsiadlowski et al. 
(2004) suggest that special evolutionary circumstances are required to lead to an explosion of 



detected gamma-ray flux of 2 — 3 x 10~^ photons cm~^ gives 20-30 photons cm~^ in 10® seconds, 
which combined with the ACIS/HMRA effective area gives total detected photons in the range 5 x 10^ — 10* 
in the remnant as a whole. 
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the hypernova type and the accompanying gamma-ray burst, presumably involving a specific 
type of binary interaction to give the required rotation rate and mass loss, van Putten 
(2004) identifies the small branching ratio with the probability of the central compact object 
remaining at the explosion center and accreting material to become a high-mass black hole. 
This argument would suggest that Cas A, which produced what is presumably a neutron 
star recoiling with velocity in the plane of the sky of 330 km s^^ (Thorstensen, Fesen, & van 
den Bergh 2001), could not have produced a gamma-ray burst. Aside from these statistical 
arguments, the recoiling neutron star should naturally be expected to destabilize the jets, 
causing rapid baryon loading and quenching any nascent gamma-ray burst. 

We have verified that the explosion center determined from the positions and velocities 
of the X-ray knots cataloged by Delaney et al. (2004) is consistent with that determined 
from optical observations by Thorstensen, Fesen, & van den Bergh (2001), although the 
uncertainties in the X-ray position are considerably larger than those in the optical work. In 
this case the central compact object recoil is at about 75 degrees to the jet axis. According 
to Wang, Lai, & Han (2005), the recoil component perpendicular to the jet axis requires 
either a relatively slow rotation rate upon explosion, rather unlikely in our view since such 
a case appears unlikely to generate jets, or the presence of a binary companion. A binary 
companion is also favored to produce the required degree of pre supernova mass loss (Young 
et al. 2005), and so this is the scenario we prefer. 

Additional circumstantial evidence against Cas A hosting a gamma-ray burst lies in the 
apparent absence of stellar wind circumstcllar media in most gamma-ray burst afterglows 
modeled to date (Chevaher & Li 2000; Panaitescu & Kumar 2002), with Piro et al. (2005) 
identifying one of the few (GRB 011121) consistent with a stellar wind density profile. Of 
course Cas A shows clear evidence of expanding into a remnant stellar wind, consistent 
with only a minority of gamma-ray bursts. The Cas A explosion also seems to have been 
extraordinarily dim, possibly due to being shrouded in a dense stellar wind at the time of 
explosion, though Young et al. (2005) comment that the dust component of the wind, if the 
Cas A progenitor evolved as a Wolf-Rayet star, is unlikely to have significantly changed the 
interstellar extinction from values found today. 

Mazzah et al. (2005) report the observation of the Type Ic supernova SN 2003jd that 
appears to be an intermediate case between normal Type Ic SNe and gamma-ray burst host 
SNe. No gamma-ray burst is noted, but the supernova is clearly asymmetrical, exhibiting 
double peaked [O I] lines. These are interpreted as coming from equatorial regions of the 
progenitors, while emission from the polar regions is dominated by jets. More recently, 
Folatelli et al. (2005) and Tominaga et al. (2005) report similar conclusions for SN 2005bf, 
another energetic Type Ib/c supernova. We suggest that the explosion of Cas A may have 
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been a similar event. It is clearly asymmetrical, with a little more mass in equatorial than 
in polar regions (taking the estimates of Table 5 in Laming & Hwang 2003, at face value), 
presumably the result of polar jets forcing the overlying stellar envelope to lower latitudes. 
However the degree of asymmetry produced appears to be at the low end of that predicted 
by models, and the jets themselves are clearly less energetic than in either the simulations 
of Khokhlov et al. (1999) or in observations of gamma-ray bursts (Frail et al. 2001). 

5. Conclusions 

We have studied the NE polar region of the Cas A SNR in considerable detail. In 
contrast, the counter jet region appears to be expanding into a more complex circumstellar 
medium, possibly a cloud of higher density, which limits the conclusions wc may draw. Fits 
to knots in the counter jet are completely consistent with those in the NE jet stem. However 
the high ionization age knots at the tip of the NE jet, which turn out to be of crucial 
importance for the analysis in this paper have no counterparts in the counter jet. It is also 
worth commenting that the three jet tip knots studied in detail here are simply not visible on 
images from the earlier 50 ks Chandra observations of Cas A. The deep 1 Ms VLP observation 
was absolutely necessary to see these at all. Analyses of the knots reveals that the "jet" -like 
morphology really is due to an explosive jet, and not arising as a result of interaction with a 
cavity or lower density region in the circumstellar medium. While Cas A exhibits many of the 
properties suggested for gamma-ray burst hosts; extensive presupernova mass loss and jets 
presumably requiring nonnegligible rotation of the progenitor, the jet itself is inferred to be 
significantly less energetic than generally accepted for gamma-ray bursts, and so is unlikely 
to have actually been a gamma-ray burst itself. Any pair or Poynting flux dominated jet is 
hkely to have become baryon loaded if the jet was destabihzed by the recoil of the neutron 
star, giving rise to the ejecta dominated structure that we see today. Further, the jet is 
also significantly less energetic than in the jet-induced explosion model of Khokhlov et al. 
(1999), suggesting that the supernova giving rise to the Cas A remnant was not a jct-induccd 
explosion. This is supported by Laming & Hwang (2003) who find evidence of asymmetry in 
the explosion at the low end of or below the range generally modeled or invoked to explain 
polarization observations of Type Ibc or Type II supernovae. Consequently we believe we 
are seeing jets produced as a by-product of an explosion that proceeded by other means, 
e.g. a more usual neutrino generated event, and that evidence of convective overturn in the 
ejecta should be interpreted in terms of these types of models (Kifonidis et al. 2000, 2003, 
e.g.) rather than the rotating jet powered models (Akiyama et al. 2003, e.g). 

This work only scratches the surface of the analysis to be done and the physics to be 
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extracted from the Cas A VLP dataset. In subsequent papers we intend to analyze in detail 
the ejecta knots, particular emphasis on the Fe rich regions. Using methods developed in 
Laming & Hwang (2003) and Hwang & Laming (2003), together with observed knot veloci- 
ties, both from proper motions in the plane of the sky and from Doppler shifts, to reconstruct 
a three dimensional distribution of ejecta knots. With this in hand, we then envisage study- 
ing various explosion models to gain insight into the mechanisms and instabilities involved 
by comparison with our catalogue of knots. 

This work was supported by grants from the CXO GO and the NASA LTSA programs. 
JML was also supported by basic research funds of the Office of Naval Research. 



A. Jet Models in BLASPHEMER 

In Laming & Hwang (2003) and Hwang & Laming (2003) we used core-envelope solutions 
given by Truelove & McKee (1999) for trajectories of the forward and reverse shocks, where 
the ejecta density distribution is taken to be represented by a uniform density core with an 
outer power law envelope p cx r~", where n > 5. Here we outline our implementation of the 
Truelove & McKee (1999) solutions for power law ejecta density distributions, p oc r"" with 
n < 3, which have no constant density core, expanding into a circumstellar medium with 
p oc r~^. For reasons given in the main text, we expect these to be better representations of 
the jet regions, taking s = 2 as in Laming & Hwang (2003). We use the same system of units 
as in Laming & Hwang (2003)^. Other parameters are; the lead factor /ed = 1.1+0.4/ (4 — s), 
the ratio of forward reverse postshock pressures, 0ed = 0.343 (1 - s/3)°■^^ and the expansion 
velocity of the outermost ejecta = 2-^ (5 — n) / (3 — n). 

The blast wave trajectory for n < 3 ejecta during the ejecta dominated phase is given 

by 

Rh 



t = 



ejI'ED 




ED r-,(3-s)/2 



-2/(3-n) 

(Al) 



We define a transition time to Sedov-Taylor behavior as follows. The blast wave radius at the 
transition is estimated from = Rb/t — 2^/^rI^~^^^'^ / (5 — s) with t — Rb/vejlED yielding 
-Rconn = (i^ej^ED (5 — s) /2yJX)^^'^^ ^\ The transition time is then calculated from equation Al 
with i^conn, where ^ = ^/(S — s) (10 — 3s) /Stt. For t < tconn, Rb is calculated from equation 



^The numerical factors in equations Al and A2 in Laming & Hwang (2003) for the s = ambient medium 
are in error; they should be 473.6 and 3.43, instead of 423 and 3.07, respectively. 
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Al, Rr — Rb/lED, 




n-3 



ED r>(3— ■s)/2 



3-S V lEDfn^b 



(3-5)/2 



and Vr — Rr/t — Vb/lsD- For t > tc 
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t 



(A3) 



where the reverse shock velocity is held constant, similarly to the core propagation phase in 
Laming & Hwang (2003). In this approximation, if R^/t < 0.75vr, the reverse shocked ejecta 
flows backwards towards the remnant center. In this case we put the expansion velocity of 
reverse shocked ejecta equal to zero, hence 



Rr — 0.75Vrt 



Rn 



^ED^conn 



— Vr In 



(A4) 



giving t = tconnexp (rconn/ZED^conn " 0.75) and so 



Rr = 0.75t'rtconn exp 



-^conn 
^ED^conn'J^r' 



0.75 , 



(A5) 



and then Vr — 0.75Rr/t. While a httle ad hoc, this gives good agreement with the s — 
models given in Truelove & McKee (1999). We emphasize that the jet ejecta in which we 
are most interested undergo reverse shock passage with t < tconn, and so the reverse shock 
velocity is given accurately by equation Al and those in the text immediately following. The 
later behavior of Vf, and Rb is only necessary for the adiabatic expansion of these shocked 
knots. 

The remaining parameter required is the separation between the forward shock and the 
contact discontinuity. In a jet, this can be significantly smaller than the value obtained in 
a spherically symmetrical explosion, because plasma entering the forward shock in the jet 
region is shocked obliquely, and its postshock flow is no longer radial. We set this distance 
equal to zero in our simulations, so that the contact discontinuity is at the forward shock. 

We have also updated the atomic data for dielectronic recombination of K-shell, and 
L-shell ions. Recombination from H- to He-like and from He- to Li-like are taken from 
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Dasgupta & Whitney (2004). The successive isoelectronic sequences Li-, Be-, B-, C-, N-, 0-, 
and F-hke are taken from Colgan, Pindzola, & Badnell (2004) (see also Colgan, Pindzola, & 
Badncll 2005), Colgan et al. (2003), Altun ct al. (2004) (see also Altun et al. 2005), Zatsarinny 
et al. (2004a), Mitnik & Badnell (2004), Zatsarinny et al. (2003) (see also Zatsarinny et al. 
2005), and Gu (2003) respectively. Additionally dielelectronic recombination from Ne- to 
Na-like and from Na- to Mg-like are taken from Zatsarinny et al. (2004b) and Gu (2004). 



B. The Kompaneets Approximation 



We give a treatment of the Kompaneets approximation for Gas A, following Bisnovatyi- 
Kogan & Silich (1995), in the case of a uniform explosion into a nonuniform circumstellar 
density distribution. Wc then comment on the modifications necessary to treat asymmetric 
explosions into symmetric circumstellar media. The fundamental assumption leading to the 
Kompaneets approximation is that for an explosion into asymmetric media, the pressure 
behind the blast wave is constant with position on the shock surface. Defining the shock 
surface hy f {r,9,t) — in conditions of azimuthal symmetry, then 



dt dt drdt dddt dt ' ^ 



(Bl) 



Writing I V/l = ^J{^f/^rf + {^f/T^ef and \v\ = ^(7 + 1) P,ft/2p = ^J {dr/dtf + (rde/dtf 
we get. 



dr 
di 



2p 



1+1 ^ 
r do 



(B2) 



We assume a density profile p = Po exp (O/W) (ro/r)*, where po and Tq are the values of p and 
r at some fiducial point, taken here to be the head of the jet, and put y — t-\J (7 + 1) Pgh/^pQ 
so that 



— — 



dyj 



1 



(91nr 

de 



exp {-e/w) . 



(B3) 



At this point wc remark that a postshock pressure variation of P = PQeic^ {—9 /W) and 
a uniform ambient density would be mathematically identical. However such a pressure 
variation needs to be sustained somehow, since otherwise with a high postshock sound speed 
the pressure behind the blast wave should equihbrate. We have already seen that in the jet 
region, shocked circumstellar plasma encounters the blast wave obliquely, and once shocked 
flows non-radially away from the jet tip. Consequently one might expect the dynamics of 
the forward shock in the jet region to remain ejecta dominated well past the time that the 
forward shock elsewhere in the remnant has evolved to the Sedov-Taylor phase. The cooling 
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time of shocked ejecta at the jet tip is also sufficiently short to render the sound speed at 
the jet tip slow, and hence impede the pressure equilibration among the shocked ejecta. 
Therefore we will assume that the expanding jet "remembers" the anisotropy of explosion 
energy well into the remnant phase and calculate blast wave profiles from the corresponding 
modification of the Kompaneets approximation. 

Wc define ^ = {ro/rY^'^ dr/dy = dro/dy, evaluated at ^ = but also assumed indepen- 
dent of 9 below, rearrange and integrate to find 



lnr = ± J VC^exp {e/W)-l d9 + ln J ^dy 



(B4) 



and hence 



dlnr 



-de 



Putting ^ = cosh x and J ^dy = Tq we can integrate to find 



\_ 2W 



(B5) 



2W I 'J 

VC^ exp (e/W) - 1 + arctan (^^C^ exp (^/W^) - l) + 



arctan 



(Ve^)B6) 



where the -ve sign has been taken. If the normal to the shock front makes an angle a to the 
radial direction, then dlnr/dO — tana and tan^a = ^^exp (O/W) — 1. 

Lazzati & Begelman (2005) give a preferred angular dependence of dE/dQ oc 9^"^ outside 
of a central region of constant energy. This is derived from considerations of the gamma-ray 
burst jet-cocoon structure at jet breakout from the stellar envelope, and also reproduces the 
so-called Frail correlation. In this case our Kompaneets approximation becomes 



Tq / dr 
\dy 



dlnr 
89 



9^ ■ 



(B7) 



Following the steps above we derive 



In^ 



2e 



9 + 



9^ 



^2 ni 



9 + J9^ 



W2 
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2^2 ni 
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n W + W^1- ^{R8) 



Simphfied expressions in the limit 1/^^— >Oor^~l and 9 »W are r — Tq exp {—^9'^/2W) 
and r = ro {29/W)^^^ ex.-p {-9^ /2W) respectively. Zhang, Woosley, & Heger (2004) prefer 
n — 3, and Graziani, Lamb, & Donaghy (2005) prefer n = 8. These cases require numerical 
integration of the corresponding version of equation B7. 
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Table 1. NEI Models for Jet Tip Knots with O Continuum 



Region 


Counts 


Region Size 


X^dof 


Nh 


kT 


net 


Si 


Fe 


N tip 


3134 


4.2 X 2.2 


78.5, 1.11 


1.2 


0.60 


1.6cl3 


2200 


640 












(0.55-0.66) 


(>1.4el3) 


(1170-4000) 


(250-1600) 


M tip 


9626 


9.5 X 2.7 


187.9, 1.36 


1.3 


0.73 


9.6el2 


2.3 


0.44 












(0.70-0.75) 


(>4.6el2) 


(1.6-2.6) 


(0.26-0.59) 


S tip 


2895 


r=1.5 


92.2, 1.40 


1.38 


0.60 


9.4cl2 


3.1 


4.2 










(1.30-1.46) 


(0.58-0.62) 


(> 6.5el2) 


(1.5-672) 


(2.1-673) 
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Table 2. Line Strengths in Jet Knots 



Region 


Counts 


Region Size 


Fe K centroid 


Fe K EQW 


Ca He a EQW 


Si Ho/Ly ratio 


kT 




North Filament 


jet3 


19953 


5.7 X 4.9 


6.632 


4.5 (3.6-5.3) 


0.99 (0.83-1.2) 


2.7 (2.0-3.5) 


0.74,3.4 


j4 


38989 


9.6 X 2.1 


6.635 


3.7 (3.1-4.3) 


0.76 (0.71-0.83) 


3.3 (3.1-3.6) 


[0.17], 2.4 


j5e 


19120 


4.4 X 4.2 


6.690 


2.4 (1.3-3.3) 


1.0 (0.81-1.3) 


4.2 (3.2-4.8) 


0.45, 3.7 


j5a 


19279 


7.5 X 3.5 


6.643 


2.6 (2.1-3.2) 


0.88 (0.68-0.82) 


2.2 (1.9-2.4) 


1.0, 11 


j5c 


32972 


7.2 X 3.2 


6.722 


2.5 (2.4-2.6) 


0.53 (0.45-0.62) 


5.7 (4.9-6.7) 


0.41, 2.1 


j5b 


27203 


8.4 X 2.1 


[6.661, 0.05] 


0.9 (0.07-1.7) 


0.19 (0.13-0.30) 


3.9 (3.0-4.7) 


0.58, 1.8 


Middle Filament 


j8 


23291 


7.6 X 4.0 


6.592 


2.8 (2.0-3.6) 


1.2 (1.1-1.4) 


6.3 (5.5-7.3) 


0.46, 3.7 


jio 


37045 


9.1 X 2.4 


6.601 


2.3 (1.5-3.1) 


1.0 (0.91-1.2) 


6.0 (5.6-6.7) 


0.55,2.6 


jll 


25017 


7.0 X 2.6 


6.666 


2.0 (1.4-2.7) 


1.1 (0.93-1.2) 


3.5 (3.0-4.0) 


0.53, 3.7 


jlla 


29189 


9.2 X 3.8 


6.626 


2.1 (1.6-2.7) 


0.64 (0.55-0.72) 


3.3 (3.0-3.8) 


0.61, 2.9 


South Filament 


jl3a 


39392 


4.6 X 2.2 


6.605 


1.7 (1.1-2.1) 


0.76 (0.68-0.83) 


4.2 (4.1-4.5) 


1.8, 80 


jl3b 


20171 


5.3 X 1.6 


[6.610, 0.1] 


1.3 (0.3-1.7) 


0.43 (0.36-0.51) 


6.6 (6.0-7.7) 




jl8 


74142 


11.2 X 2.9 


6.645 


0.58 (0.44-0.73) 


0.10 (0.08-0.12) 


4.0 (3.8-4.3) 


1.7, [150] 


jl9 


62144 


8.0 X 2.8 


6.686 


1.2 (1.0-1.4) 


0.17 (0.15-0.20) 


3.1 (3.0-3.3) 


1.7, [150] 


j20 


44787 


6.0 X 2.3 


6.676 


1.7 (1.3-2.1) 


0.22 (0.18-0.28) 


2.8 (2.7-3.0) 


1.6, [100] 


j21 


60072 


9.7 X 2.6 


6.648 


1.3 (0.95-1.5) 


0.11 (0.09-0.14) 


2.7 (2.5-2.8) 


1.6, [150] 



Counter-Jet 



cjet 60518 



6.597 4.0 (3.0-5.1) 0.37 (0.30-0.45) 4.4 (3.4-5.3) 0.39, 1.86 
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Table 3. pshock Spectral Fits for Jet Knots with O continuum 



Region gain factor x^/dof Nh kT n^t Si Fe 



North Filament — pshock 



jets 


1.001 


254.4, 


l.SO 


1.57 


2.32 


2.1ell 


1.2 


0.68 










(1.51-1.64) 


(2.09-2.56) 


(1.8-2.6ell) 


(0.96-1.5) 


(0.56-0.90) 


j4 


0.999 


338.4, 


1.54 


1.29 


2.00 


2.4ell 


1.1 


0.58 










(1.28-1.32) 


(1.89-2.13) 


(2.2-2. 7ell) 


(1.0-1.2) 


(0.52-0.67) 


j5b 


0.997 


224.3, 


1.32 


1.18 


1.26 


5.6ell 


0.56 


0.23 










(1.13-1.21) 


(1.15-1.48) 


(3.4-9.0ell) 


(0.48-0.61) 


(0.19-0.26) 


j5e 


0.999 


194.8, 


1.11 


1.20 


1.74 


2.5cll 


0.96 


0.52 










(1.14-1.21) 


(1.59-1.94) 


(2.1-2.9cll) 


(0.84-1.0) 


(0.40-0.55) 


j5c 


0.999 


428.2, 


2.25 


1.18 


1.40 


2.4ell 


1.1 


0.50 


j5a 


0.999 


347.7, 


1.61 


1.17 


2.79 


2.1ell 


0.68 


0.28 










(1.14-1.20) 


(2.51-2.96) 


(1.9-2.5ell) 


(0.58-0.72) 


(0.24-0.30) 


North Filament— NEI 


jets 


1.001 


350.6, 


1.79 


1.65 


2.15 


l.lell 


2.1 


1.0 










(1.59-1.70) 


(1.93-2.35) 


(9.4cl0-1.2cll) 


(1.2-18) 


(0.75-38) 


j4 


1.000 


670.8, 


3.05 


1.25 


1.85 


l.lell 


1.1 


0.51 


j5b 


0.997 


273.2, 


1.61 


0.86 


1.07 


5.1ell 


0.45 


0.91 










(0.77-0.96) 


(1.03-1.15) 


(3.5-6.4ell) 


(0.38-0.49) 


(0.07-0.12) 


j5e 


0.999 


275.9, 


1.57 


0.81 


1.15 


3.8ell 


0.54 


0.13 










(0.75-0.85) 


(1.13-1.21) 


2.7-4.5ell) 


(0.46-0.58) 


(0.10-0.14) 


j5c 


1.001 


693.4, 


3.65 


1.03 


1.24 


1.2ell 


0.73 


0.24 


j5a 


1.003 


589.0. 


2.73 


1.04 


3.55 


e.OelO 


0.49 


0.18 


Middle Filament — pshock 


j8 


0.998 


240.0, 1.26 


1.49 


1.96 


1.5ell 


1.5 


0.58 










(1.45-1.54) 


(1.91-2.17) 


(1.4-1.6ell) 


(1.3-1.9) 


(0.50-0.70) 


jio 


0.999 


285.7, 1.43 


1.28 


1.62 


1.8ell 


1.8 


0.69 
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Table 3 — Continued 







/do f 


Nh 


kT 


n t 


Si 


Fe 








(1.26-1.30) 


(1.53-1.79) 


(1.6-2.3ell) 


(1.5-1.9) 


(0.52-0.71) 


jll 


0.998 


254.8, 1.32 


1.22 


1.83 


2.7ell 


1.5 


0.50 








(1.16-1.26) 


(1.68-2.03) 


(2.3-3.4ell) 


(1.2-1.8) 


(0.41-0.61) 


jlla 


1.000 


270.3, 1.26 


1.26 


1.86 


2.5ell 


0.73 


0.26 








(1.21-1.28) 


(1.76-1.98) 


(2.2-3.1ell) 


(0.66-0.76) 


/ (\ C\C\ (\ C\C\\ 

(0.20-0.28) 








Middle Filament — 


-NEI 






j8 


0.998 


315.3, 1.66 


1.33 


1.85 


7.8el0 


1.1 


0.33 








(1.28-1.39) 


(1.58-2.00) 


(7.0el0-1.0ell) 


(0.90-1.2) 


(0.22-0.36) 


jio 


0.998 


421.4, 2.11 


1.20 


0.99 


4.2ell 


2.0 


0.57 


jll 


0.998 


348.4, 1.81 


0.86 


1.20 


3.9ell 


0.89 


0.15 








(0.80-0.92) 


(1.16-1.26) 


(3.2-4.6ell) 


(0.78-1.0) 


(0.13-0.20) 


jll 


0.999 


375.3, 1.75 


0.79 


1.26 


3.4ell 


0.51 


0.066 








(0.78-0.81) 


(1.24-1.37) 


(3.2-3.9ell) 


(0.49-0.52) 


(0.057-0.074) 


South Filament — pshock 


jl3a 


1.006 


358.4, 1.66 


1.39 


2.12 


l.Scll 


1.7 


0.26 








(1.35-1.43) 


(2.05-2.30) 


(1.3-1.6ell) 


(1.6-1.9) 


(0.23-0.31) 


jl3b 


1.009 


229.6, 1.26 


1.50 


1.96 


l.Oell 


1.2 


0.24 








(1.47-1.55) 


(1.79-2.13) 


(9.2el0-l.lell) 


(1.1-1.3) 


(0.20-0.26) 


jl8 


1.008 


493.6, 1.61 


1.47 


1.95 


1.8ell 


0.28 


0.11 








(1.46-1.49) 


(1.85-1.98) 


(1.7-1.9cll) 


(0.27-0.29) 


(0.10-10.) 


jl9 


1.008 


518.9, 1.87 


1.37 


2.32 


l.Scll 


0.50 


0.19 








(1.35-1.39) 


(2.22-2.42) 


(1.7-1.9ell) 


(0.47-0.52) 


(0.17-0.20) 


j20 


1.008 


450.8, 1.91 


1.21 


2.27 


2.1ell 


0.74 


0.18 








(1.19-1.25) 


(2.11-2.40) 


(1.9-2.4ell) 


(0.68-0.78) 


(0.15-0.20) 


j21 


1.008 


421.9, 1.60 


1.20 


1.68 


4.3ell 


0.59 


0.15 








(1.17-1.22) 


(1.61-1.76) 


(3.9-4.9ell) 


(0.57-0.63) 


(0.13-0.16) 



South Filament— NEI 
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Table 3 — Continued 



Region gain factor /dof Nh kT ngt Si Fe 



jl3a 


1.005 


469.5, 2.16 


1.21 


2.08 


7.3el0 


1.1 


0.15 


jl3b 


1.004 


308.3, 1.69 


1.36 


1.85 


5.8el0 


0.88 


0.14 








(1.31-1.40) 


(1.73-2.04) 


(5.3-6.7el0) 


(0.79-0.92) 


(0.11-0.15) 


jl8 


1.007 


731.9, 2.39 


1.27 


2.02 


7.3el0 


0.23 


0.075 


jl9 


1.008 


807.6, 2.92 


1.20 


2.64 


6.6el0 


0.39 


0.13 


j20 


1.008 


771.9, 3.27 


1.02 


2.87 


6.6el0 


0.56 


0.10 


j21 


1.008 


567.7, 2.16 


0.72 


1.58 


2.0ell 


0.37 


0.037 



Counter- Jet — pshock 



cjet 1.000 315.7, 1.36 2.06 1.26 4.1ell 1.22 1.2 

(2.04-2.11) (1.24-1.29) (3.8-4.5ell) (1.17-1.43) (0.97-1.43) 



Counter-Jet— NEI 



cjet 1.000 427.6, 1.57 1.95 1.08 2.9ell 0.79 0.63 

(1.93-2.03) (1.05-1.09) (2.6-3.9ell) (0.74-1.04) (0.56-1.00) 
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Fig. 2. — Spectra of the outermost jet knots (north at upper left, middle at upper right, 
south at bottom). 
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Fig. 3. — The Fe K equivalent width plotted for knots in the north (black), middle (red), 
and south (green) jet filaments as a function of the distance of the center of the extraction 
region relative to the position of the Cas A point source near the center of the remnant. The 
correlation of Ca EQW compared to Fe EQW is shown in the panel to the right. 
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8.50 9.50 10.50 11.50 12.50 13.50 

log(ionization age s/cm3) 

Fig. 4. — The loci on the electron temperature - ionization age plane of fitted knots in 
the NE jet, color coded such that knots from the middle filament are red, knots from the 
south filament are green, and those from the north filament are black. Also plotted are 
curves derived from models or circumstellar cavities for ejecta profiles corresponding to a 
uniform density core and a power law outer envelope, p oc r~", with = 7, 9, and 11. The 
discrepancy between observations and models suggests that CSM cavity models cannot be 
responsible for the jet morphology. 
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8.50 9.50 10.50 11.50 12.50 13.50 

log(ionization age s/cm3) 

Fig. 5. — The same fitted jet knots as in Figure 4, compared with jet models specified to 
match the fitted abundances in the knots at the jet tips, as sohd fines with the same color 
coding. Small differences in the onset of thermal instability are visible, due to the different 
abundance sets. The black dashed line shows the result of assuming half the plasma in the 
N filament to be H dominated. This reduces the electron temperature by over a factor of 
two, but the onset of thermal instability is less sudden. Jet models appear to match the 
data points much better. 
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Fig. 6. — Charge state distributions of Fe (top) and Si (bottom) corresponding to ejecta 
encountering the reverse shock at different times following explosion. The thick solid line 
indicates Fe (or Si) going through the reverse shock f .6 year after explosion, which gives the 
best agreement with the fitted electron temperature. The dashed curves with sucessively 
finer and finer dashes are for ejecta undergoing reverse shock passage 1.45, 1.4, 1.35, 1.315, 
and 1.3 years after explosion, with electron temperatures 7.84 x 10^, 7.04 x 10^, 6.17 x 10^, 
5.51 X 10^ and 5.22 x lO^K respectively. Although these ejecta cool to lower temperatures 
than indicated by the fits, the charge state distributions are in better agreement with the 
observations. The narrow solid lines indicate the collisional ionization equilibrium charge 
state distributions for a temperature of 8.6 x 10^ K, and active indicative of the actual 
ionization balance detected in the jet knots. 
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